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Abstract. We present results from a 20ksec RXTE ob- 
servation of the black hole candidate Cyg X-l. We apply 
self-consistent accretion disk corona models to these hard 
state data and show that Comptonization in a spherical 
corona irradiated by soft photons from an exterior cold 
disk is able to successfully model the spectrum. We also 
present the power spectrum, the coherence function, and 
the time lags for lightcurves from four energy bands. By 
modeling the high-resolution lightcurves with stochastic 
linear state space models, we show that the rapid hard 
state variability of Cyg X-l can be explained with a single 
timescale. 



1. Introduction 

The X-ray emission of galactic black hole candidates 
(BHCs) displays different "states" cha racterized by dis- 
tinct temporal and spectral properties ( Esin, McClintock 
|fc Narayan, 1997; van der Klis, 1995| ). Except for occa- 
sional transitions to the soft state, Cyg X-l spends most 
of its time in the hard state (Fig. |l|). It is therefore the 
prime candidate for evaluating accretion models of this 
fundamental black hole state. 

Since 1996 the Rossi X-ray Timing Explorer (RXTE) 
provides data from 3 to 200 keV with millisecond time- 
resolution. In this paper we derive a wide variety of spec- 
tral and temporal hard state properties by applying phys- 
ical accretion disk corona models (ADCs) (Sect. ^) as well 
as new timing methods (Sect. [| and ||) to a RXTE hard 
state observation of Cyg X-l. Whereas ADCs are very suc- 



cessful in explaining the observed spectrum (Dove et al 
1997), a unified model for the hard state accretion also 
has to take into account the time dependence of the emis- 
sion. Some of the associated constraints that are provided 
by the power spectral density, the coherence function, the 
time lags, and the lightcurve modeling with stochastic pro- 
cesses, are discussed in sections Ej and For further dis- 



cussion see 



et al. (1 



Nowak et al. (1998j 1999|) and |Pottschmidt 
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Fig. 1. RXTE- ASM count-rates in the 1.5 to 12keV band 
for Cyg X-l from the beginning of the ASM measurements 
until 1997 December. The arrow indicates the time of our 
pointed RXTE observation. The 1996 soft state can be 
clearly identified by the increased count-rate from 1996 
May through August. 

2. Data Extraction 

RXTE observed Cyg X-l for a total integrated source time 
of slightly more than 20ksec, starting 1996 October 23, 
briefly after its transition to the hard state (Fig. |l|). 

We analyzed data from the Proportional Counter Ar- 
ray (PCA) as well as from the High Energy X-ray Timing 
Experiment (HEXTE). An overview of our data screen- 
ing procedures and of the RXTE instruments is given by 
Dove et al. (1998) ) and Wilms et al. (this volume). For the 
spectral analysis we used 3 to 30keV PCA data and 20 
to 200 keV HEXTE data. The temporal analysis is based 
on PCA lightcurves with 1.95 ms resolution from four en- 
ergy bands of approximately equal count rates (0(700- 
1000)ctss _1 per band). The energy bands are defined in 
Fig. ||, the total length of the lightcurves was ~18ksec. 
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sphere+disk geometry where the seed photons come from 
an accretion disk external to the inner spherical corona, 
allows to fit the broad-band spectrum of Cyg X-l with 
one self-consistent model. 

The best-fit spere+disk spectrum for our RXTE obser- 
vation is shown in Fig. ^|. The deviations between the data 
and model above lOkeV are consistent with the current 
calibration uncertainty of the PCA, while the deviation 
in the iron band between 5 and lOkeV shows the phys- 
ical simplifications in our model. This does not pose a 
severe problem to the model since the deviations repre- 
sent less than 1% of the total flux modeled. The deviation 
might be due to an underestimation of the covering factor 
for reflection, e.g., by a possible overlap between the disk 
and the corona due to the transition from the soft to the 
hard state two weeks before the observation (cf. Fig. [j]). 
The presence of accretion disk flares (due to disk heating) 
would also increase the reprocessing observed. Finally, a 
slight increase in the Fe abundan ce will also produc e the 
stronger line seen in the data. See Dove et al. (1998 ) for a 
discussion of these points, as well as for a presentation of 
spectral fits using phenomenological "standard" models. 



Fig. 2. Fit of the sphere+disk model to the data, indicat- 
ing good overall agreement between the ADC model and 
the observations (x 2 c d = 1-55, assuming a systematic error 
of 2% for the PCA data). The corona has a temperature of 
kT= (65.7±3.3)keV and an optical depth of r = 2.1 ±0.1. 
The seed photons were assumed to come from an accre- 
tion disk external to the corona, with a r -3 / 4 temperature 
profile and an inner accretion disk temperature of 150 eV. 

More details on the PCA timing data of this observation 



are given by Nowak et al. (1999) 



3. Spectral Modeling: Self-Consistent Thermal 
Accretion Disk Corona Models 

The BHC hard state spectrum, a power-law with T ~ 1.65, 
modified by an exponential cutoff with Sfoid ~ 150 keV is 
naturally explained by accretion disk corona (ADC) mod- 
els, i.e., by Comptonization of low-energy seed photons 
from a cold, geometrically thin accretion disk in a hot, 
semi-relativistic corona (Titarchuk, 1994). As we explain 
in more detail in an accompanying paper (Wilms et al., 
this volume; see also Dove, Wilms & Bcgelman, 1997 



Dove et al., 1997), we have developed an ADC model 



in which the radiation field, the temperature and opac- 
ity of the corona, and the reprocessing of coronal radia- 
tion in the accretion disk are calculated self-consistently. 
The commonly assumed slab geometry where the accre- 
tion disk is sandwiched between two coronae is not self- 
consistent, since the coronal temperatures and opacities 
are not within the the physically allowed range of com- 
binations of these parameters. On the other hand, the 



4. Timing: Fourier Methods 

Power Spectral Density 

We calculated power spectral densities (PSDs) by divid- 
ing the data into contiguous segments of uniform length 
and computing the Fourier transform Sj(f) of each data 
segment. The discrete PSD is then given by (15^ | 2 ) where 
the brackets indicate an average over data segments and 
logarithmic frequency bins (Fig. see also Nowak et al. J 



1999) 



The PSDs show characteristic hard state behavior (cf. 
Belloni fe Hasinger, 199C| ). For 10- 3 -100Hz the rms vari- 
ability in all bands is about 30%. The PSDs are relatively 
flat for frequencies from 0.02-0.2 Hz. For higher frequen- 
cies, they are roughly described by a steepening, broken 
power-law with a break at ^2 Hz, and typical slopes of 
— 1.0, and —1.5, respectively. Note that while there are no 
strong energy dependences, the high frequency slope does 
"harden" for higher energies. 

4-2. Coherence Function 

The coherence function, 7 2 (/), is a frequency-dependent 
measure of the degree of linear correlation between two 
concurrent time series (Vaughan & Nowak, 1997; Nowak 
]pt al., 199S| ) . Specifically, it gives the fraction of the vari- 
ability of one time series that can be predicted from 
the other. The coherence is based on the Fourier trans- 
forms S(f) and H(f), calculated for the soft and the hard 
lightcurve, respectively: 



7 2 (/) 
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Fig. 3. Energy resolv ed PSDs normalized according to 
Miyamoto et al. (1992| ). PSDs for the higher energy bands 
have been multiplied by powers of 10 for clarity. The Pois- 



Fig. 5. Time lag as a function of frequency between the 
lowest energy band and the three higher on es. The error 



son co unting- noise has been subtracted taking the instru- 
mental deadtime into account (Zhang et al., 1995). 



bars a re calculated according to eq. (16) of Nowak et al 
(1999|) . 
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PSDs. The coherence is remarkably constant and unity 
over the broad frequency range from 0.02 to 10 Hz (Fig. |4|) 
- in agreement with previous results and contrary to the- 



oretical expectations (Vaughan & Nowak, 1997; Cui et al 



1997 ). This behavior of indicates that the Comptoniz 



ing medium must be static on these timescales, despite 
the various possibilities for non- uniform coronae (Nowak 



et al., 199£ 
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Fig. 4. Coherence between the lowest energy band and the 
three higher ones, corrected for Poisson counting-noise. 
The noise correction and the error bars have been calcu- 



For all energy channels there is a trend for the co- 
herence to drop below ^0.02 Hz and above ^10 Hz with 
the deviation from unity coherence generally becoming 
slightly greater for increasing energy difference. The loss 
of coherence below 0.02 Hz is consistent with the viscous 
timescale for R > 50GM/c 2 , while the loss of coherence 
above 10 Hz is consistent with dynamical timescale for 
R < 50GM/c 2 (R is the distance from a 10 M Q black 
hole). For further discussion of the physical implications 
see [Nowak et al. (1998|) and |Nowak et al. (1999[). 



lated according to eq. (8) of Vaughan fc Nowak (1997 ) 



Especially for high frequencies the noise-subtracted co- 
herence is subject to additional systematic uncertainties 
(<30% from 30 to 100 Hz), which is why deviations to 
greater than unity can occur. 

where the brackets indicate an average over data segments 
and frequency bins and the star denotes complex conju- 
gation. 

We computed the coherence function between the low- 
est energy band and the three higher ones using the same 
lightcurve segments and frequency rebinning as for the 



4-3. Time Lags 

The Fourier time lag is a frequency-dependent measure of 
the time delay between two time series. Like the coherence 
it is calculated from their Fourier transforms S(f) and 
H(f), 



r(f) 



arg(S(/)*g(/)) 
2tt/ 



(2) 



t can be either positive or negative. For our sign conven- 
tion, a positive time lag indicates that the hard lightcurve 
lags the soft lightcurve. 

We computed the time lags for the same energy bands, 
lightcurve segments, and frequency bins as used for the co- 
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herence (Fig. ||). Within the frequency range that is not 
noise-dominated (0.1-30 Hz), the data are consistent with 



the so fjtest energy band always leading the harder energy 
bands. The lags increase significantly with the energy dif- 
ference between the bands and become larger for smaller 



process. This is consistent with our results for model- 
ing E XOSAT observations of Cyg X-l ( pottschmidt et aT 



1998) and in analogy to the latter work we derive a re- 
laxation timescale r = (0.13 ± 0.03) s for the RXTE 
lightcurves. In contrast to the usually applied shot noise 



frequencies, a result previously found by Miyamoto et al. model fits in the frequency domain (Lochner, Swank & 



(1992 |)j The time lags range from ~0.002 s to ~0.05 s Szymkowiak, 199l| ), the LSSMs of first order do not re- 



and can be very roughly described with a power law 
T if) K 7 > superposed by characteristic breaks at ~0.5 
3, and 10 Hz. 

The simplest expectation in the framework of Comp- 
tonization models is that time delays should be due to dif- 
ferent diffusion times through the corona and thus should 
be nearly independent of Fourier frequency ( Miyamoto & 
Kitamoto, 1989). Furthermore, the overall energy depen- 



dence of the time lags is also not compatible with current 
Comptonization models ( Nowak et al., 1999 ) - including 
the sphcrc+disk model. It is therefore nec essary to con- 
sider more complicated accretion models (Nowak et al., 
1998[)~[rhe smallest observed time lags might indicate that 
the Comptonizing medium is rather small, < 30 GM/c 2 . 
Consequently, if the time delays are due to disturbances 
propagating through the corona, the propagation speeds 
are extremly slow (O(l~10% c)). These values are incon- 
sistent with current ADAF models. 



quire pre-defined shot forms and are able to describe the 
rapid hard state variability of Cyg X-l with one temporal 
parameter. 

The derived relaxator can in principle be used for scal- 
ing time dependent ADC models, by fitting LSSMs to the- 
oretical "ADC lightcurves" . First estimates indicate much 
larger coronal radii than implied by the shortest time lags 
(~ 600 GM/c 2 ). For a better understanding of spectral 
und temporal hard state properties in one unified accre- 
tion scenario, we plan to apply the models and methods 
presented here to further RXTE observations of Cyg X-l 
and other BHCs. 
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5. Timing: The Linear State Space Model 

An alternative method to the timing analysis in the fre- 
quency domain is to examine the X-ray variability directly 
in the time domain by modeling the observed lightcurves 
with linear state space models (LSSMs) (Konig & Tim 
mer, 199~7| ) . These models describe the dynamics of a 



lightcurve y(t) with a stochastic, autoregressive (AR) pro- 
cess x(t) of order p (Scargle, 1981): 



(3) 



where e(t) is a Gaussian white noise processes with zero 
mean and variance a 2 . The dynamical parameters, aj, are 
closely related to the temporal parameters characterizing 
the system. An AR[p=l] process, e.g., is described by a 
single relaxation timescale r = — l/ln(ai). 

LSSMs also take the observational noise into account. 
For an AR[1] process 



y(t) = c ■ x{t) + V (t) 



(4) 



where r](t) is the observational noise, again given by a 
Gaussian white noise processes with zero mean and a vari- 
ance of a 2 , and where c is a constant. 

We fitted LSSMs of different order to 842 contiguous 
lightcurve segments (0.0-14.1 keV) with a length of 16 s 
each. To evaluate the goodness of the fits, a Kolmogorov- 
Smirnov test for white noise residuals was performed. We 
find that the lightcurves can be modeled with a LSSM[1] 
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